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ABSTRACT 

This report describes the activities during the fourth six month period 
of the investigation of acoustic propagation in the atmosphere with a real- 
istic lapse temperature profile. A significant error was detected since the 
■previous semi-annual report and has been corrected in both the plane wave and 
point source solutions. This report then describes both of these problems in 
some detail along with presenting some numerical results from the model. 
Work will begin this summer on the model of propagation in an inversion. 
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1. INTRODUCTION 

A significant error was discovered in the formulation of the mathe- 
matical solution to both the plane wave problem and the point source 
problems following the previous semi-annual report [1]. The discovery of 
•this error required the reformulation of both problems and invalidated much 
of the previous work. The reformulation has been carried out and is 
described in detail in this report. Of particular interest are the values 
to be used of multi-valued functions. The incorrect use of multi-valued 
functions, in part, resulted in the error described above. Great care must 
be taken in both the formulation and the evaluation of the solution to the 
problems described below. 

The point source problem is, of course, the problem of most interest 
but the plane wave solution is most easily physically understood of the two 
and thus allows the clearest reasoning concerning the multi-valued functions 
which occur in the problems. For this reason the plane wave problem is 
discussed first and extensively. 

Parts of the work described here were carried out by graduate students. 
Mr. Greg Landheim dealt mainly with numerical evaluation of the plane wave 
problem. Numerical evaluation of the saddle point approximation to the 
point source problem was carried out by Mr. Alex Cheng. Mr. Ma Yiping is 
currently improving the basic saddle point approach with particular concern 
to the pole that occurs near the path of integrat'^on. This pole contrib- 
utes a surface wave-like behavior to the solution. In addition, Mr. James 
Brown is currently working on a method of predicting the location of the 
caustic, or shadow boundary, for the point source problem which can be 
used in conjunction with an experiment that is anticipated to be conducted 
this summer. 
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2, PLANE WAVE SOLUTION 


The problem of progagation of plane waves In a lapse temperature 
gradient is of interest for three reasons: 1) the point source solution is 
made up of the superposition (inverse Hankel transformation) of plane waves 
.over all angles, 2 ) because the plane wave solution only makes physical 
sense at real angles It is much easier to interpret and consider various 
possible values of the multi-valued functions involved, 3) the plane wave 
solution is of interest for its own right in that it models acoustic 
propagation far from the source. 

The problem to be considered is governed by the equation 

1 32p 

v2p . 0 (1) 

(z) 3t^ 

where the sound speed a|^(z) is given by 




1 + aZ 


)) 


( 2 ) 


Equation (2) implies that temperature as a function of height above the 
ground is given by 

Tm = T. + at — ~ (3) 

1 + aZ 

In addition the surface boundary condition 

- Zw = P (4) 

must be applied on z - 0 where the vertical component of the acoustic 
fluid velocity, w, is obtained from 

3w 1 3P 

3t p 3Z 

* ^ 

The solution must behave as 
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(1) w 

i(li)t --y X cos 0 +-y- Z SlO 6) 
p X e “ • 


(1) (Jj 

i(o)t - -g~ X cos 0 - -j- z sin 0) 
R(0)e “ 


( 6 ) 


as z This condition requires the solution behave like plane waves 
incident at an angle 0 in an homogeneous atmosphere when outside the region 
of strong gradients (i.e., z •♦• «). 

Taking 

i (d)t - -a- X cos 9) 

P = F(z)e - (7) 
and substituting into Equation (1) yields (where the subscript «» has been 
dropped) 


d2p 




(• 


1 + az 


- COs2 e) F = 0 


dz2 a2 1 + az + AT/T 
a second order differential equation with a turning point at 


( 8 ) 


-TP 


1 AT 

— {•—ctn'^ e - 1). 

a T 


( 9 ) 


Nayfeh [2] gives a method of obtaining an approximate solution of this 
equation valid when X = w/a a >> 1. (Note that x when used in this report 
is not the wavelength but is approximately the reciprocal of wavelength 
nondimensionalized by a.) This requirement is the same as found for the 
governing equation, Equation (1), to be valid (Reference [3]) so no new 
approximations are required. The approximate solution can be expressed as 

= — - hi (n(z)) + . - h? (n(z)) (10) 

^ g"(z) / g'(z) 

where 

3 2/3 

n(2) = (j X) ' 9 (z), (11) 
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and 


hi(5) ■ h|J] (IjS/Z) (12) 

h2(0 « (|)l/3 5 I /2 (13) 


Values of the modified one-third order Hankel functions, hj and h 2 » are tab- 
ulated and their properties discussed in Reference [4]. The functions hi 
and h2 can also be represented in terms of Airy functions as 

hi(0 = k [Ai(- 0 - i Bi(- c)] (14) 
and 

h2(0 - k* [Ai(- O + i Bi(. K)1 (15) 
where k = (12)^/^ and k* is the complex conjugant of k. The func- 
tion 


2 / sin^ e (1 + az + aT/T) - (aT/T) 

g'(z) = -a 


g(z) (1 + aZ + AT/T) 


(16) 


giving 


/ o AT ^ AT AT 

g3/Z(z) - J sin^ e (1 + aZ + — )2 - _ (1 + aZ + — ) 


1 AT 1 


2 T sin 0 


1 + (^ 

tn (- ) 

1 - ^ 


where 




sin2 6 (1 + az + AT/T) - AT/T 


(17) 


(18) 


sin^ 9 ( 1 + az + AT/T) 

This solution is rather complex and serveral important features need to be 
discussed to understand it. 

The functions hi(n) and h 2 (n) have a complicated behavior (Reference 
[4]). For real values of n both hi and h 2 yield a complex result which 
is oscillatory with an algebraic decay of the amplitude for increasing 
magnitude of the argument. The function hi can be found to represent 

downward propagation waves, incident waves, and h 2 upward propagating 


T r 


waves* reflected or refracted waves. The oscillatory behavior also occurs 
for hi when the phase of n equal to 2ir/3 and for h 2 when the phase of n 1s 
-2n/3. When the phase of n Is rr/3, hi decays exponentially and h 2 grows 
exponentially. When the phase is -it/ 3, h 2 decays exponentially and hi 
;grows. 

The function g^/^(z,B) has two branch points. The first occurs when 



and corresponds to the turning point height for a given Incident angle 6. 
At this branch point the magnitude of g3/2 is zero. A second branch point 
occurs at cos e equal to 1. Real values of cos e greater than one, of 
course, do not correspond to real angles. The magnitude of g3/2 is infinite 
at this branch point. The discussion of the branch points in terms of 
cos 9 rather than e alone will be useful in the point source problem and is 
thus used here. When cos e is real and between zero and the first branch 
point the observer, located at a height z, is the region “illuminated" by 
the incident waves. For values of cos 0 real and between the first and 
second branch points the observer is located below the turning point 
(and caustic) and is in the "noni 1 luminated" or shadow region. From 
physical considerations then, when the values of cos e under consideration 
fall between zero and the first branch point an oscillatory solution 
would be expected. When the values fall between the first and second 
branch point the receiver is in the shadow and an exponentially decaying 
solution is to be expected. 

This physical reasoning is supported by the mathematics in that for 
the values of cos e between zero and the first branch point g3/2 is real 
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and positive and a possible value of g (and therefore n) 1s also real and 
positive yielding the oscillating (propagating wave) solution. When cos e 
Is real and between the first and second branch points g^/2 is Imaginary. 
If the Is chosen to be -1 as Is generally the case when using el**^^ as 
♦‘opposed to then g^/2 is real and positive for cos e between zero and 
the first branch point and Imaginary and positive between the first and 
second branch points. The function g (and therefore n) Is then real and 
positive between zero and the first branch point and Is chosen to be complex 
with a phase of n/3 between the branch points for real cos e. Thus hi and 
hg are complex oscillating functions between zero and the first branch point 
and h 2 1s exponentially decaying while h 2 is exponentially growing between 
the branch points. 

Thus from a mathematical point of view as the observer descends from 
Infinite height toward the turning point the value of g (and n) approaches 
zero monotonically through positive real values. At the turning point 
height the value of g Is zero corresponding to the first branch point. 
Descending below the turning point the magnitude of g is a monotonically 
increasing function with a constant phase of tt/S, Since both hi and h? 
are entire functions no discontinuities occur in this process. The solution 
we are desling with, Equation (10), however includes /g^, where g^ is the 
derivative of g with respect to z, and since g and g^ are real for values 
of cos e between zero and the first branch point (turning point) and complex 
between the first and second branch points g' (and /p) are discontinuous 
at the first branch point. Thus the two terms in Equation (10) are also 
discontinuous, but the total solution must be continuous (no abrupt changes 
in either the magnitude or phase of the acoustic pressure are allowed). 
Thus different forms of solution must occur above and below the turning 
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point. It is possible, however, for the turning point to occur at a nega- 
tive value of z. This means that the turning point theoretically occurs 
below the ground surface, and in reality reflection from the ground occurs 
before refraction can turn the waves upward. Thus three formii of the 
solution occur; the reflected case where the observer is above the turning 
point, the refracted case where the observer is below the turning point, 
and the refracted case where the observer is below the turning point but 
above the ground. The first form is obtained by requiring the solution to 
satisfy the surface boundary condition, (4), and that the solution take 
the form of (6) as z •♦• ». The second two forms of the solution are obtained 
by requiring the form of the solution valid below the turning point to 
satisfy the surface boundary condition, (4) and the form valid above the 
turning point to satisfy (6) and the two forms to be continuous at the 
turning point. The process gives the values of A and B in Equation (10) as 

AT -1/2 

^ K for 0 < cos 0 < (1 


A = 


AT -1/2 

K for cos e > ( 1 + “) and z > zjp 


Ke 


TT 

l“6 


( 20 ) 


,n 

i'c 


AT 

for cos e > ( 1 + — ) 


- 1/2 


and 0 < z < zjp 


e ® - iR 


AT -1/2 

K • R for 0 < cos 6 < (1 + — ) 

T 


B = 



AT 


- 1/2 


for cos 6 > (1 Y") and z > zjp 


(21 


AT 


- 1/2 


for cos e > (1 +~) and 0 < z < ZJP 


e'6 - iR 


- 7 - 


J 


K « s1n^/^{e) 


;^l/6 ^1/2 e^(5w/12) 


( 22 ) 


t hi (n(0)) ■¥ lij> hi^ (n{0)) 

R « (23) 

T h 2 (n(0)) + i\\i h 2 " (n(0)) 


1 2 g"'(0) 

2 Po Cq 9'*(0) 


(24) 


♦ = ^2/3 g-(0) (1)2/3 (26) 

Po Co 2 

This solution differs from the forms presented earlier [1, 10] in that the 
matching at the turning point was not required in the earlier solution. 
This error produced a continuous solution but with a cusp at the turning 
point. Figures 1 to 3 present the results of evaluation of the solution 
presented above. Further discussion of this solution and the results are 
given in Appendix I, which is a paper scheduled for publication by the 
Journal of the Acoustical Society of America in September, 1984. 


3. POINT SOURCE 


The point source problem also contained the same error as the plane 
wave problem and required reformulation. As described in [3] the problem 
under consideration is governed by the wave equation, 

1 a2p 


V2P = Q 


(26) 


but including a source term 


«(r) . 

Q = q 6(z - s) ( 27 ) 

nr 


and the speed of sound profile is given by 
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1 

1 + aZ^ 


(28) 


The same profile as considered In the plane wave prable^t* Equation (2), 
The boundary condition on z * 0 Is also unchanged from the plane wave 
problem and Is given by 

- Zw « P . 

As z + a radiation condition, requiring all waves to propagate upward 
must be applied. 

The solution is obtained by first separating out the time dependence 
by defining a function l3(r,z) such that 

P(z,r,t) = ^(z,r). 

Taking the Hankel transformation (two-dimensional Fourier transform) of the 

resulting governing equation and the boundary condition yields 

d^G (1)2 1 -I- az q 

p2] G = - — 6(z - s) (29) 

dz^ a'^ 1 + az + AT/T 2 tt 


and 


6 * 


iZ 8G 

uipQ 9z 


(30) 


on z - 0. The function G($,z) is the Hankel transform of ^(r,z) and is de- 
finded by 

G(e,z) *= / G(r,z) r Jo (8r) dr (31) 

0 

The inverse Hankel transform is defined by 

G(f'.t) - / G(8,z) 3 Jo (Pr) d$ (32) 

0 

Equation (29) in the regions 0 < z < s and z > s and the boundary 


condition, Equation (30), are the same as for the plane wave case provided 
that cos e in the plane wave case is replaced by (a/o))3. Thus the same 
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solution, Equation (10), applies here. 

The plane wave problem required three different forms depending on 
the location of the receiver In physical space and the angle of Incidence 
of the wave at Infinite height. This problem also requires several forms 
but further complications occur. In the plane wave problem all real 
'incidence angles were possible, and cos 6 could vary between zero and one. 
In the point source problem the range Is more limited. Noting that 3 Is 
the equivalent of (w/a) cos e where e Is the Incidence angle of the plane 
wave at Infinity, a physical Interpretation of the mathematical solution is 
possible. Consider first waves propagating downward from the source located 
at height s. Waves that are directed downward at a sufficient angle will 
reach the ground as In the plane wave esse. Plane wave rays directed 
downward at an infinite height with incidence angles between 90® and some 
value 00, or 0 < 0 < 3 q = (w/a) cos 6© = (oi/a) + aT/T)] can pass 

through the source location and reach the ground where they are reflected 
(see Figure 4). The turning point of these waves are below the ground 
surface. The initial angle 6o characterizes the waves whose turning point 
is at the ground level, z = 0. Waves in this range of initial angle and 
their reflections from the ground constitute two groups of possible waves 
and one of the forms of the solution. 

A second group of waves which are initially directed downward at the 
source are rays with angles at Infinity in the range between e© and 
or 00 < ^ < ^z ‘ ®z * (w/®) *'[(1 + oz)/(l + Qtz + AT/T)3. The 

angle 0z characterizes the wave whose turning point is located at the ob- 
server height i4 This range of angles or of beta then characterizes waves 
that pass through the source location but whose turning point Is above the 
ground level but below the observer (see Figure 5). These waves, before 
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passing through the turning point, ^nd after passing through the turning 
contribute a second form to the solution. 

A third group of waves propagating downward at the source have angles 
at Infinity In the range between 0 ^ and 65 or < e < 0s “ (“/*) 05 ® 

(ui/a) / [(1 + as + AT/T}]. The value 65 characterizes the wave with Its 
'’turning point at the height s. It is the wave that leaves the source hori- 
zontally and Is the limiting case of waves leaving the source and propa- 
gating downward. Waves in this group have their turning point above the 
observer, see Figure 6 , and thus cannot contribute an oscillating term to 
the total solution but rather have an exponential behavior. This behavior 
Is further discussed In connection to the saddle point method of evaluation. 
This group of waves and their continuation after being refracted upward 
constitutes a third form of the solution. 

The fourth form of the solution is due to waves which leave the source 
and proptjH>: upward. These waves have 65 < 6 < 90® or 0 < 0 < 3s ^re 
not of the same form as described earlier as their amplitude must be correct 
for waves leaving the source rather than waves reflected or refracted upward. 
These waves constitute a fourth form of the solution. 

All of the waves of the four types described above can be seen in the 
ray diagram of the point source Figure 7, but a number of additional types 
of waves need to be included in the total solution which do not appear in 
the ray diagram. The fifth form of the solution has 65 > 9 > 0® or 3s 3 
< (i)/a, these are waves with their turning point above the source and do not 
appear in the ray diagram but contribute to the full solution. In this case 
g( 0 fS) is a complex number with phase of n/3 and the solution has an exponen- 
tial behavior. 

The actual forms of the solution are obtained by requiring (10) to 
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satisfy the surface boundary condition, (30), at 2 ■ 0, at z * s, the source 
height, the solution must satisfy the condition that 


lim [G(3 ,z)3 « lim Cg(p,z)3 

2-^S. Z-fS+ 


which requires continuity of pressure and the condition that 


36 

lirn M 
z-^s 


11 m [-^3 

Z.S, 


q_ 

2ti 


(33) 


(34) 


which Is obtained by integrating Equation (29) from z = s-eto 2 = s+ e 
and taking the limit of e 0. At the turning point height the condition 

11m [6(3, z)3 = lim [ 6 ( 3 , z)3 (35) 

Z+Z]'p_ Z-+-ZTP4. 

is required if the turning point is above ground surface, z *= 0. The 
resulting solution for z > s is given by 

6 = K h 2 (n(3»z)) [hi(n(3,s)) + Rq h£ (n(3,s))3 (36) 

for (ij/a > ^ 3s^ 3o’^ 3 ^0 

6 = K h 2 (n(3,z)) [hi(n(3,s)) + R\ h 2 (n(3,s))3 (37) 

for ui/a > 3z > 3s > 3 > 3o > 0 

•JL 

G = K hg (n(S,r)) Rj [hi(n(B.s)) + II 2 (n(S,s))] (38) 

for ii)/a > 3z > 3 > 3s > 3o > 0 and 

.JL .“f' 

G = K[e ^ hj(n(B.z)) + Rq hj (n(B.z))] Ri [^i (n(S,s)) 

+ Ro h 2 (n(6,s))] (39) 

for o)/a > 3 > $2 > 3s > 3o > 0* When s > z the forms are 

6 - K h 2 (n(3,s)) [hi (n(3»z) + Rq h 2 (n(3»z))3 (40) 

for u)/a > Bs > :3z > 3o > 3 > 0 

6 “ K h 2 (n(3,s)) [hi(n{3»z)) + Ri hi (n(3»z))3 (41) 

for o)/a- > 3s > 3z ^ 3 > 3o > 0 
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(42) 


G * K h 2 (n($.s)) [h^ (n(e.z) + Rq hg {n(3»z))] 

for u)/a > 0s > e > 02 ^ ^0 ^ ® 

G = K Ce’^ hj (n(e,s)) + U2 (n(0.s))] e'^Rj [hj (n(e»z)) 

+ Ro ^2 (n(6,z))] (43) 

for o)/a > 0 > 0s > 02 > 00 > 0» where 



T hi (ti(0,O) + 1 > hi' (ri(0,O)) 

Ro r 

T h2 (n(0,O)) +i4^ h2" (ri(0»O)) 


and T and 41 are defined in (24) and (25), These eight forms account for 
the five physical situations described above since some forms are repeated. 
The solutions include both direct waves, waves which have not been reflected 
or refracted and the reflected or refracted waves. For 0 < 02 the inter- 
pretation of the many forms are aided by recalling that hi (n(0»z)) repre- 
sents downward traveling waves and h2(n(3,z)) represents upward traveling 
waves. For 0 > $2 the interpretation is more complicated as the terms 
have an exponential behavior (hi growing and h2 decaying with increasing 
z). The entire discussion above is based on real values of beta, but beta 
must be interpretable as a complex number and at this point the boundaries 
between the various forms of the solution off the real axis must be clar- 
ified. 

A square root occurs in the function g3/2(0,z) and the branch line for 
this square root must be chosen so that 1 = - i as mentioned in the 
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plane wave section (see also Clemmow [5]). If the branch line for the 
square root is chosen to be along the line where the phase of the argument 
of the square root is tt/ 2 then the desired result is obtained. This results 
in branch lines for the function g^/2(e,z) extending downward from two 
branch points (3 = Pz or e = w/a) in the real 3 > 0 half plane and upward 
in the real ( 3 ) < 0 half plane. These branch lines for g^/^(0,O), g^/2(3.z) 
and g^/2(3»s) then connect directly to the boundaries on the real axis and 
form natural boundaries when real ( 3 ) > 0 and imaginary ( 3 ) < 0 where the 
solution is discontinuous. The argument in the modified Hankel function 
involves g(3.z) etc., rather than g3/2(3,z) and a branch line is necessary 
for the two-thirds power also. Here, as was pointed out in the plane wave 
solution (i)2/3 must be e‘'’^/3. If the branch line for the two-thirds 
power is chosen to lie along the line where the phase of is tt 
then the desired result is obtained. This branch line in the 3 space is a 
curve from the first branch point of g^/^(3.z), 3 = (w/a) ✓ [(1 + azT/(l + 
az + AT/T)], to the second branch point m/a and located in the imaginary 
( 3 ) > 0 region on the real ( 3 ) > 0 half plane (see Figure 8). These lines, 
one for each g(3»0) g(3.z) and g(3»s) constitute the remaining parts of the 
boundaries of the various forms of the solution. Thus it appears that their 
are eight different forms of the solution each discontinuous from the 
others. This can be shown to not be the case. 


On crossing each of these branch lines the phase of the function g 
jumps from - Ztt/S to 2n/3 and the phase of gz discontinuously increases 
by 27 t/ 3, Reference [4] gives some useful results for this type of behavior: 



.2tt .n .7? 

ho(n (e’3 hj (n) + hg (n)) 



(47) 

and 

,2tt 

hj (n e^^) = - h 2 (n) 

(48) 
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using these results it is easy to show that 

Rj(r,) • Ro (n e ’3 ) (49) 

which in turn demonstrates that equations (36) and (37), and (40) and (41) 
are continuous at the branch line for g(3,0) ■ (g^/^(3»0))2/^. Similarly, 
, using Equations (47) and (48) it is possible to show that (37) and (38), 
and (41) and (42) are continuous across the g(3,z) branch line and (38) 
and (39), and (42) and (43) are continuous across the g(3»s) branch line. 
This exercise in showing continuity also points out another way of writing 
the solution, equations (36) through (43), These modified forms of the 
solution are particularly useful in that the first and second term in each 
is continuous with the respective term in the forms of the solution 
bounding it, and the first term always represents the direct wave and the 
second the reflected or refracted wave. For 2 > s the solution can be 
written 

G = K |^h 2 (n(3,z)) [hi (n(3»s) + Rq (n(3»0)) h 2 {n(3,s))]j (50) 

in region A of 3 space (see Figure 9) 


. 2tt 

G = K |h 2 (n(3,z)) [hj (n(3.s)) + (n(3,0) e’^ j (n(3,s))]| 

in region B 

f . ZTT 

h 2 (n(3,z)) [hi (n(3,s) e^^ ) 

.2il .2tl 

+ Rq ^n(3,0) e^2 ) h 2 (n(3,s) e^^ )]| 

r -ijL 

G = K h 2 (n(3.z) e’^ )] hi (n(3.s) e^^ ) 

^ ..Zjl 

+ Rq (n(3,0) e’3 j (n(3,s) e’^ )] 


(51) 


in region C 


(52) 


(53) 


and in region D and for 2 < s 


- 15 - 


» 


I 


G * K |^h2 (n(0.s)) [hi (n(0,z)) + Rq (n(3»0)) hg (n(6.z))3j (54) 

In region E 2^ 

G . K |^h 2 (n(8.s)) [hi (n(B.z)) + Ro (n(B,0) hj (n(B.2))]j(55) 


In region F 


2n 


G = K hg (n(3,s)) C hi(n(3.2) e'*’ ) 


2ti 


+ Rq (n(3,0) e’^ ) hg (n(e»z) e'*^ )] 


2tt 

I3 


(56) 


In region G 


2tt 


2tt 


G = K j h 2 (n(3,s) ) [h^ (n(3,z) 

^ .2tt .2tt ”) 

+ Ro (n(B.O) e’3 ) h 2 (n(B.z) e’^)]j 


(57) 


and In region H. This then completes the formal solution to the problem 
with the solution given In the form of an Integral, the Inverse Hankel 
transform, Equation (32), where the function G(0,z) is given by Equations 
(50) through (57). This Is a very complex Integral and approximations 
are necessary to proceed toward Its numerical evaluation. 


4. APPROXIMATE INTEGRATION OF THE POINT SOURCE SOLUTION 

The formal solution discussed In Section 3 Is extremely complex and 
there Is very little hope of carrying out the required Integration exactly. 
Thus an approximate integration of the solution is required. Approximate 
solutions of propagation problems of this type are discussed by Keller 
[6], Because of the physical insight obtained and the success of the 
approach In the homogeneous atmosphere propagation problem, the saddle 
point approach, or the ray approximation, as termed by Keller, has been 
attempted. 

The saddle point method is a method of approximately integrating 
expressions of the form 
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J F(B) eEf(B) dB 


(B8) 


where 5 is a large parameter. If we interpret this integral as an integral 
in complex p space, the path of integration can be modified according to 
the rules of contour integration in the complex plane. What would appear 
to be desired is then to find a point (or points) in the complex p plane 
where the integrand of (58) is a maximum and to distort the path of inte- 
gration to pass through this point in such a way that the exponential term 
in (58) decreases rapidly as one moves along the modified path of integra- 
tion away from the point of maximum integrand. In the complex plane, 
however, a regular function does not have a maximum, but rather a point 
where 

3f 

T7 = 0 (59) 
B=Bsp 

is a saddle point of the function, a point where the function surface 
looks like a mountain pass, col, or saddle. This point is found to give 
the major contribution to the integral when the path is modified to give 
a rapidly decreasing integrand as one moves away from it. Approximately 
integrating (58) then can be shown to give 


/-i" f(6sp) 

/-Af"(Bsp) 

where o is the angle at which the path of integration crosses through the 
saddle point. This method is described in more detail and with more 
rigor in most texts on advanced mathematical methods, for example [7,8]. 
This method is not a fool-proof one however, as the modified path of 
integration may intersect branch lines, pass around poles, or lead into 
other difficulties. Keller [6], however, asserts that the major part of 
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the solution is obtained. Also, if no saddle point occurs then the inte- 
gral can be approximated as zero. Using this rather crude approach the 
Integrals can be approximated relatively easily. 

To use the saddle point approach the formal solution (32) and (50) 
through (57) must be put in the form of (58). The change in the limits of 
’integration is easily achieved by using the identity 



Jo(re) = j[ + Ho<2)(rs)] 

(61) 

and noting that 

pb) 

(62) 

and that 

G{B,z) = G(- B.2) 

(63) 

yieldi ng 

00 



Q =1 J G(a,z) (ar) da 

mOO 

(64) 


The exponential forms are obtained by then assuming the arguments of all 
the Hankel functions involved are sufficiently large to allow replacement 
by their asymptotic forms, given in [4] and [9]. The resulting expressions 
are of the form (58) with the integrand of each integral consisting of two 
terms, one will be seen to represent the direct waves, and the second the 
reflected and/or refracted waves. This becomes clear when the derivative 
with respect to B is taken of the argument of the exponential terms to find 
the saddle points. It is found that 

3(x g3/2{3,z)) 

= - F(3,z) (65) 

33 

and the resulting expressions are exactly the expressions obtained for 
the rays of the point source problem, section 2, with the substitution 
of cos 6 = (a/w) 3 so that the rays are expressed in terms of a. Thus, the 
finding of the saddle points, which could be numerically quite difficult, 
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has a clear physical interpretation in terms of rays cf the point source. 
For an observer located at a point in physical, (r,z), space which is 
"illuminated" by the point source, the saddle points of the approximate 
Integrals giving the solution correspond to the 0 values of the two trays 
passing through that point. Figure 10. 

At any point in the "illuminated" region of physical space the solution 
is then approximated, as described above, by two of the terms listed 
below. For direct waves, with z > s, Region 1 and 2 of Figure 10 


^ ^-1(Sspr+n> 

^ 



(66) 

where 


i>l(3) = - X g^/2(3,s) + X g^/2(s,z) - j 

(67) 

and 0SP is the root of 


r + (0 ) * 0 

^1 sp' 

(68) 


For direct waves with 0 < z < s. Regions 3 and 4 of Figure 10 

-i(6sprH2) 


gd = 


/ 2jf K 2 G 


✓ - ( 3 r 

^2 SP 


where 


♦ 2 ( 6 ) = X 9^''^ ( 6 ,s) - X 9 ^^^ (B,z) 


(69) 


(70) 


and 3sP is the root of 


r * <®sp) ' ° 


For reflected waves. Regions 1,3 and 5 of Figure 10 

Rp(Bc„) 


Gr = 


^2 

/■: TTp 7 

^3 SP 


(71) 


(72) 
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where 


^3 * ^ 9^^^ (3»s) + ^ (0»2) “ 2 X (3sp» 0) ■ J (73) 

and 3sp is the root of 

r •♦■ ♦a (PSp) ® 0 (74) 

For refracted waves which have not contacted the caustic, Region 4, 5 and 6 


of Figure 10 


where 


Kg e’’ f®SP''+*4) 

-: (3SP) 


/ + 4> Cb T 

^4 SP 


(3) = X g^/^ (3sp»s) + X (Ssp.z) +y 


,3/2 


11 


and 3sp is the root of 


r + ♦; (3SP) = 0 
4 


(75) 


(76) 


(77) 


For refracted waves that have contacted the caustic, Regions 2, 4, and 6 of 


Figure 10 


where 


Gr 


/ 2tt K£ e 


-i (3 sp^+<I>5) 


✓ - 6" (3 J 

^5 SP 


(78) 


<t»5 ^ <l>4 (79) 

and 3sp is also a root (and smaller of the two roots for region 6) of 
(77). The derivative of g3/2 is defined by (65) and is given in Reference 
[ 10 ]. 


K2 (Bsp) " 


and 


2BsP 1 


1 


1 


1 


12i Xtt / irr v'g^TF^^TzT /q^TBsp^ 9^'^^(BsP*2) 9^'^^(6sP*5) 


3 1/3 

T + i|» i (7 X) /g(3,0) 

R2(3) * ;a- ,r/3 

T - ij; i (7 X) /glM) 


(80) 


(81) 


Thus the problem has resolved itself into a relatively simple expres- 
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sion, but two difficulties remain. First, the expressions given above are 
not continuous as the receiver is moved from region to region in physical 
space, and secondly, no expression has been obtained that is valid in the 
shadow, region 7, Figure 10. The first difficulty is easily resolved by 
an ad hoc approach. This discontinuity results from the fact that on the 
boundary between regions the arguments of one of the modified Hankel func- 
tions is zero and it has been assumed that all of the arguments are large 
enough to allow use of the asymptotic forms. The actual integrals are 
continuous at these boundaries but the asymptotic forms are not. Thus by 
replacing the asymptotic forms of the modified Hankel functions in (66), 
(69), (72), (75) and (78) by the modified Hankel functions themselves one 
obtains in place of (66) 

% = — hg (n(Bsp»2) hj (n(8,s)) (82) 

/-!?"■> 

in place of (69) 


Gq = 


’^'^3 .-i(Bspr- f) 


hi (n(0sp*^)^ ^2 (’^(^SP*^^^ 


93 


(<i'2) 


in place of (72) 

^ / 2ti K3 Rq^^sp^ -i(3sp^-7) 


Gr = 


^ / V 

- J^(*3) 


hg (n(3sp»z)) hg (n(35p,s)) (84) 


in place of (75) 

/ 2tt K 3 Ri(3^p; -i(05pr-7r) 


Gd = ^ e 

— 

A rrz (<i>4) 


93 


hg (n(35p»z)) hg (n(35p*s)) (85) 
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and in place of (78) 

Gr “ — ^ J) (n{35p,z) h£ (n(8sp»®)) W 


a continuous result is obtained. Here 

q 1 


Kg = 


2Psp 


(87) 


24 i i/g^JJ^z) /g^TPs^) 

The second problem, an expression valid in region 7, the shadow, can 
also be circumvented. Examination of (85) and (86), the expressions valid 
in region 6, shows that these expressions become infinite at the caustic, 
where 




t5 

8e2 


= 0 


( 88 ) 


This result is easily understood by realizing that the refracted rays are 
given' by r = F(3,s) + F(e,z), Equation (77) expressed in F, and that the 
caustic is defined by the largest radius that can be obtained for fixed 
observer height, z, and the source height, s, or where 3r/3p = 0, yielding 
the above expression. This singularity, in the mathematical sense, results 
from expanding the argument of the exponential terms in the saddle point 
method in a Taylor series but retaining only second order terms, which are 
zero in this case. In the physical sense it arises from the ray tube going 
to zero and conservation of energy. A result valid near the caustic has 
been obtained by Sachs and Silbiger [11]. Using their approach, which 
essentially involves expanding the exponential terms about the ray tangent 
to the caustic rather than the saddle point, yields the expression 


Gr = 2^ K 3 (eJ Ri(Bc) ^ (( ir) 


♦ (Be) 

h 2 (n(Bc.z)) *>2 (n(6c.s)) 


♦ (Be) 


(89) 
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where Be the root of 


(0c) * 0 

(90) 

which is valid near the 

caustic 

and on both the 

"illumlnsted" and the 

shadow side of it. Here 

Ar 

* r - re 

(91) 

where the caustic radius 


" (0c) 

(92) 


The solution has not been carried out into the shadow region beyond the 
region where (89) is valid at present. The method to do so is clear how- 
ever. For observer locations in the shadow, saddle points exist and near 
the caustic have a real part approximately equal to the value of B at the 
caustic. Since the saddle points are no longer on the real B axis as they 
were in the "illuminated" region, they are no longer physically interpret- 
able as in terms of angles, but give rise to the "complex rays" discussed 
in the underwater acoustics literature [11, 12], 

5. RESULTS 

A limited number of cases have been calculated and plotted using the 
procedure described above. This work has been mainly carried out by Alex 
Cheng and Yiping Ma. The cases presented below are all for values of AT/T = 
0.025, a = 2.5 m“l, s = 3 m, the same conditions as the ray diagram of 
Figure 10, and mainly for a frequency of about 1600 Hz (o) = 10,000 rad/sec). 
Variation of aT/T, a, and s appear to have little effect on the general 
nature of the solution other than to move the interference pattern and 
shadow boundary* in physical space. Figures 11 to 15 are plots of sound 
pressure level (with an arbitrary reference) versus elevation at radii 
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of 20, 40 and 80 m. Figure 11, r •« 20 m, is sufficiently close to the 
source that most waves directed toward the surface reach It and are re- 
flected, only waves leaving the source at very shallow angles are refracted 
upward before reaching the ground, see Figure 10. The regions of Figure 10 
are indicated on Figure 11 for reference. For most of the vertical range 
'presented an Interference pattern with a mean sound level of about 32 IB 
occurs. The oscillations in this interference pattern slowly Increase in 
amplitude with height except for heights near the source height, s 3 m, 

where low sound levels occur. This peculiar behavior near the source 

height and at small radii Is characteristic of the solution. Careful 
examination of a ray diagram with rays emitted by the source at equal 
angular if^tervals show that the rays are spread very far apart in vertical 
distance in this region of space. If the equal acoustic energy is assumed 
to be present between rays, a location with a large distance between 

rays should be a region of low acoustic pressure and sound level. This 

appears to be occurring in the region near the source height and to account 
for the low sound levels seen there. The slowly increasing amplitude of the 
oscMlations is due to the increasing strength of the reflected waves present 
as the height incre«ises. 

Figure 12 is a similar plot fejt at r *= 40 m. Above the source two 
maxima occur in the interference pattern above the dip occurring near the 
source height. The mean sound level in this interference Is about 26 dB, 
6 dB less than the mean on the previous plot, as might be expected. Near 
the source height occurs the dip due to the ray spreading and below this is 
the maximum that occurs just inside the shadow boundary. At a height of 
about 1.6 m a discontinuity occurs in the sound level. This results from 
the approximation near the caustic which is not required to be continuous 
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with the other forms of the solution, the program, however, attempts to 
make the solution continuous if possible. 

Figure 13 is again a plot of height versus sound level. Note that the 
height scale does not continue to the ground at z « 0. Again the regions 
of Figure 10 are shown on the plot. Here discontinuities in slope occur 
at about 4,2 and 22 m. The first is at the boundary between the dircfct 
and the refracted waves and the second is between the approximation near the 
caustic and the other forms of the solution as discussed above. Again the 
riiean sound level has decreased about 6 dB to near 20 dB with the doubling 
of the distance from the source as compared to the previous figure. 

Figures 14 and 15 are at a radius of 40 m but at higher frequencies, 
Figure 14 at w » 12000 rad/sec (f *= 1910 Hz) and Figure 15 at w « 15000 
rad/sec {f = 2390 Hz). Although plotted at different scales the figures 
show very similar behavior to Figure 12. 

Figure 16 is a plot of sound level versus horizontal distance from the 
source. Again the regions of Figure 10 are shown along with the empirical 
model of Wiener and Keast [13]. The Wiener and Keast model suggests that 
a simple source model (corrected for atmospheric absorption which is not 
present here) yields the sound level between the source and the shadow 
boundary. Beyond the shadow boundary an empirically obtained excess atten- 
uation is applied to the simple source model. This empirical model, based 
on the experimental data they collected shows good agreement with the local 
mean sound level in the interference region and good agreement with the slope 
of the decay of the sound level in the shadow. Note that a simple source 
cannot be expected to produce the interference pattern but only the local 
mean sound level. The Wiener and Keast model would be in better agreement 
with the present model in the shadow if the excess attenuation was applied 
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after the radial location where the sound level falls below the simple 
source value rather than at the actual shadow boundary which always has a 
sound level greater than the simple source does at the boundary. Figure 
17 is a comparison of the excess attenuation seen In this case as compared 
with the model of Wiener and Keast. Note that if Wiener and Keast did not 
correctly determine the shadow location the •'constant" difference between 
their model and the present model would be accounted for. Thus the best 
experimental data currently available which is summari'’ed in the Wiener 
and Keast model appears to be in good agreement with the present model. 

6. STATUS 

A model of acoustic propagation In a lapse temperature gradient above 
a finite impedance ground surface has been developed and yields reasonable 
results when compared to the empirical model of Wiener and Keast, This 
model was numerically implemented by Alex Cheng, The model contains numer- 
ous approximations, however, and the experiments planned to be conducted in 
the Fall of 1984 by NASA Langley will be helpful in assessing its validity. 
Work is still underway on some improvements to the model which will allow 
evaluation further into the shadow region. Also the pole which gives rise 
to a surface wave-like term is being Investigated. Methods for locating the 
pole and determining the value of resulting residue have been developed but 
the question of when the resulting term is to be included is yet to be 
resolved. It is hoped that this question can be answered by the end of the 
summer. The question of the pole and the surface wave-like term Is being 
investigated by Ma Yiping. 

The development of a set of programs to be run on a hand-held, HP-41, 
or a VAX computer to locate the shadow boundary has been completed and 
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supplied to NASA Langley. These programs were written to support the 
experiment mentioned above which in turn will assist the modeling effort. 
It is anticipated that this work and some additional wcrk on the acoustic 
rays in the model temperature field will form the basis of an additional 
paper to be prepared. This work has been mainly carried out by James 
Brown, 

Following the completion of this report the author will begin the 
formulation of the inversion condition problem. Although this problem is 
somewhat more complex than the lapse condition problem experience with 
the simpler problem is anticipated to be extremely useful in developing 
that model . 

Finally, a considerable amount of thought has been given to the 
physical interpretation of a surface wave and a model that may give some 
insight into that question has been developed. The mathematical solution 
to this physical model is hoped to be investigated during the upcoming 
year. 
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Fig. 2. Height versus sound pressure level for: (a)y = 10, AT/T« - 0.02, 
6 = lo'», o = 2 (b) T = 10. AT/T_ ' 0.02, 6 = 2°. a = 2 ni'l . 
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Figure 4. Ray diagram for a plane wave passing through the point source 
location and being reflected from the ground. 
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Figure 5. Ray diagram for a plane wave passing through the point source 
location and being refracted before reaching the ground but 
below the receiver. 



Figure 6. Ray diagram for a plane wave passing through the point source 
location and being refracted before reaching the ground and 
above the receiver. 





Figure 7. Ray diagram for a point source showing the plane waves indicated 
in Figures 4 through 6 and others. 
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Figure 9. Regions in complex 6 where the for.s of the solution given 

Figure ( 50 ) through (B7) hold. 
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Figure 10. A ray diagram showing the seven regions in physical space with 
a = 2.5 m-1, aT/T = 0.025, and s = 3 m. 
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The behavior of a plane wave reflecting from a flnite impedance surface with a realistic 
atmospheric temperature profile is investigated. An approximate solution has been implemented 
on a digital computer and this solution is presented graphically for a number of cases at low i 
incidence angles, j 
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INTRODUCTION 


‘mMmiiN 


The atmosphere is not isothermal under normal condi* 
tions, yet most acoustic propagation models make that as- 
sumption. Temperature inhomogeneity can occur as ran- 
dom fluctuations resulting from convection cells and 
turbulence or as slow diurnal variations in the atmosphere 
immediately adjacent to the ground surface. It is the latter 
case that this paper addresses. 

During daylight hours when sufficient insolation is 
present, the atmosphere is typically warmer near the ground 
surface than far a^ve it, a lapse condition. Although these 
regions of strong temperature gradients extend, at most, 
only a few meters above the ground, most acoustic receivers 
are also in this range of altitude. Ground mounted micro- 
phones, which are sometimes used to avoid ground reflec- 
tions, are unfortunately located so as to receive the maxi- 
mum effects of refraction due to temperature gradients. 


I. TEMPERATURE PROFILE 


Any study of the effects of thermal gradients on near- 
earth sound propagation requires the selection of a tempera- 
ture profile which is both realistic and mathematically trac- 
table. Geiger’ gives some experimental results which 
indicate that the variation of temperature with altitude be- 
haves as the logarithm of the height above the ground. Ten- 
ekes and Lumley^ give an argument based on turbulence 
theory that leads to the same conclusion. 

Logarithmic profiles are difficult to deal with math- 
ematically, Most previous work (e.g., Pekeris’) on the prob- 
lem of acoustic propagation in a stratified atmosphere have 
considered the profile 


( 1 ) 


t4(f 

which yields^pie exact solution jpm the resulting differen- 
tial equation. The vertical temperature profile used here 


r. r. VI -for/ 


( 2 ) 


is similar to Eq. ( 1 ) in its behavior near the ground (z as 0) in 
that it allows a steep gradient, but, more realistically, asymp- 
totically approaches a constant value high above the ground. 
Using data obtained by Willshire^ and by Butterworth’ and 
fitting Eq. (2) to this data, typical values for or ranged from 
1.2 to 5.9 m- ' and for A T/T^ from 0.0032 to 0.0342 for 
lapse conditions. Similar values of a and magnitude of 


AT/T^ were obtained for inversion conditions. In other 
eases, particularly where a lapse condition is being replaced 
by an inversion or vice versa, the fit can be very poor. 


II. ANALYSIS 


Consideration of the problem of propagation of plane 
waves in a thermally inhomogeneous atmosphere yields con- 
siderable insight into the more complex problem of propaga- 
tion of waves from a point source. The problem is also inter- 
esting in its own right as a model of the propagation of waves 
from a point source when the source is very far from the 
refracting region. 

The problem considered is governed by the equation 


e^(z) dt^ 

with the sound speed c{z) given by 
AT 


(3) 


At the surface {z^^) the boundary condition 


(4) 




(5) 


— Zw = p, 

where w is the'vertical component of the acoustic fluid veloc- 
ity and Z the acoustical impedance of the surface, is applied. 
The solution is assumed to behave as 


P = exp i^t — -—x cos ^ -I- — z sin ^ j j 

■f /? ((? )exp|/^<a? ^ X cos 0 ^ z sin 0 )] 1.1 


II • ) 

as z -*- 00 . This condition implies the solution behaves as 
plane waves in a homogeneous atmosphere outside the re- 
gion of strong gradients. The subscript oo , indicating evalua- [ 
tion as z-^ 00 , is omitted below^to simplify notation. The 
subscript zero will be used to&dicate evaluation at the 
ground surface. '' ' ( 

Taking _ ^ 

P « F(z)cxp{ I )x cos ^ ] ) (7) 

and substituting into 1^. (3) yields S vJiLz>^ 


dtp c* V 1 -i-az-H 


( 8 ) 


cos*^V=0, 

+ AT/T ) 

II II 

. a second-order differential equation with a turning point at 


0002 


TELLTALE 


AUTHOR’S NAME 


0002 


* ^ir,>m{\/a)l{/iT/T)coim^e^ll (9) 

Nayfeh* gives a method of obtaining an approximate aolu* 
tion of this equation which is valid when y — a>/(oc)> 1. This 
requirement is the same as that found for Eq. (3) to be valid, 
so no new approximations are required, The approximate 
solution can be expressed as 

where 

(II) 


% o<cos^<(i-h4r/r)-''*, 

A, cos0>(i+iir/r)“‘'* 

and z>Zrm, 

, ( 12 ) 

and 0<z<z„, 

'A>R, 0<cos$<(li- AT 

jg ~ t*-/*) 

c»e>{l+AT/T)-''^ 
and i>x„. (13) 

- 0 ^. cosfi>(i+dr/r)-« 

and 0<z<Zrr, 


A — >Jrr/3 sin'^‘ 

r/»,(77(0)) + itph\ivm 

^ ’■A2(T7(0}) + /#i(»7{0))’ 

r=A- ±- 

2/^oCog'(0)’ 


g^'^(z) = sin^ ^ ^ 1 -f az -f 


_IAT i j 
2 r sin 0 V 1 - 

(18) 

. sin‘0(\ +az + AT/T)-AT/Ty'^ 

^ V sin^tf(l -l-az-f dr/r) J ’ 

(19) 


(20) 

and 



(21) 


Values of the modified one-third order Hankel func- 
tions, hi and hj, have been tabulated and their properties 
discussed.^ For the numerical analysis the functions hi and 
hj have been expressed in terms of Airy functions, Ai and Bi, 
where 



( 22 )' 

and 

Aa(^)-^*lAi(-f)-H/Bi(-^)],’ (23) 

where k k (12)'^^ e " and A: * is the complex conjugate of 

k. 

The asymptotic expressions for the modified Hankel 
functions are available,^ and if the parameter y is sufficiently ; 
large the usual result for the reflection of plane waves from a 
plane surface in a homogeneous atmosphere is obtained. 

The form and behavior of the solufhn depends upon 
whether the waves which are initially tr«>veling downward, 
toward the ground, are reflected upward from the ground or 
if refraction turns the waves upward before the ground is 
reached. If cos ^<[(1 the turning point is be- 

low the ground surface and reflection from the ground sur- 
face occurs. In the opposite case the turning point is above . 
the ground surface and the wave is refracted upward before 
the ground is reached. , 

Based on meteorological data whered T/Tranged up to 
0.0342, rays with incidence angles less than about 10* could ; 
be refracted upwards before reflection from the ground 
could occi^. Although these angles are quite small, they are , 
within the 'range of common **look angles” for aircrift noise ' 
testing and refractive effects can be expected to play a signifi- 
cant role in this case. 

If the waves are refracted into upward propagation at 
the turning point, the solution depends upon whether the i 
receiver is above or below the turning point. Above the turn- ! 
ing point the solution describes an interference pattern simi- 
lar to that in a uniform atmosphere. Below the turning point : 
exponential growth or decay occurs for large values of 77 . 

The first term in Eq. (10) decays as the receiver moves 
downward, away from the turning point, if V ~ 1 is taken to ! 
equal — 1 * and, similarly, the acoustic field below the turning . 
point is reduced as 0 goes to zero and the turning point ' 
moves upward, if ln( — 1) = —iir is chosen. The second i 
term in ^. ( 10 ) also decays with decreasing 0, but grows as z 
decreases from the turning point value. This behavior occurs 
since the second term represents the reflection from the 
ground surface of the exponentially decreasing first term. 

111. RESULTS 

Numerical evaluation of the solution for the propaga- 
tion of plane waves in a temperature gradient has been per- | 
formed on a digital computer. Figure 1 is a typical example ’ 
and is intended as the referent as parameters are varied in 
other drawings. Surface impedance values given by Ches- ! 
sell’s* relationship for grass have been used, with specific 
flow resistance chosen as 300 cgs units. A value of y 10 
corresponds to a frequency of approximately 1000 Hz with 
a » 2m~'. Although this value of and the additional value > 
discussed below are not extremely large they are given as ! 
examples and experiments will be required to demonstrate | 
the applicability of the theory of these values. Figures 1 to 3 ' 
show plots of height above the ground versus sound pressure ; 
level in decibels. The reference level is assumed equal to the ' 
amplitude of the incoming wave at an infinite height. Thus 
the plots show the difference between the sound pressure . 
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level of the incident wave in an isothermal atmosphere which 
is neither reflected nor refracted and the levels resulting 
from the refraction/reflection process in a nonisothermal 
atmosphere. I 

For the conditions used to calculate Fig. 1(a), f ■■ 10, 
AT /T'm 0.2, 6 tm 5*, and a»l m~ ' the turning point and 
caustic is located at a height of about 0.8 m. Note that the 
maximum sound pressure level occurs somewhat above this. 
Above the turning point an interference pattern occurs with 
a slowly decreasing sound level in succeuively higher maxi- 
ma. Since the ratio of amplitudes of the incident and refract- 
ed waves is near unity, values of the sound levels above ap- 
proxiiM^W 6^ must be due to increased amplitude of the 
waves WM^h^alues occurring far above the ground. The 
similar amplitude of the two waves also leads to the strong 
cancellation observed in the interference pattern. Below the 
turning point a rapid decay occurs as the shadow region is 
entered. 

Figure 1(b) is for the same conditions but with an isoth- 
ermal atmosphere, J T/T — 0. In this case the reflected wave 
is somewhat weaker than the incident wave and this is appar- 
ent from the levels reached in the interference maxima and 
minima. Comparison of Fig. 1(a) and (b) clearly points out 
the increased sound levels and variable spacing of the inter- 
ference pattern in the thermally inhomogeneous case. 

Figure 2(a) shows a similar case but with the incidence 
angle increased to 10*. In this case the turning point would be 
below the ground surface and reflection at the ground oc- 
curs. The reflected wave is weaker than the incident wave in 
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FIG. 2. Height venui sound pressure level for; (a) y > 10, .d T/T^ > 0.02, 
^-|(r,a-2 m“';(b)y»10,dr/r. -0.02,^-2’,a-2 m"'. 
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FIG. 1. Height versus sound pressure level for: (a) y- 10,dr/r. -0.02, FIG. 3. Height versus sound pteaiuiv level for: (a) y- I0,dr/r. -0.01, 
^«5*,o-2m'-';(b)y- lO.dr/r. -0,tf-r,a-2 ro‘'. [L tf-y,flr-2 ri“*;|b)y- I0,dr/r, -0.02,tf-5*,o- 1 m"'. 
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(his case and this is indicated by the levels occurring in the 
interference pattern. Figure 2(b) again shows a similar case 
but with ^ 2*. In this case the turning point is at a height of 

about 7,7 m above the ground and the rapid decay below the 
turning point is more apparent than in Fig. 1(a). 

In Fig. 3(a) jr/r has been decreased to 0.01 yieldinga 
turning point at about 0. 1 S m. Comparison to Fig. 1(a) shows 
that this change has reduced the envelope of the sound level 
maxima, and that the interference minima are at a higher 
sound level, both due to a slightly weaker reflected wave 
than in Fig. 1(a), but stronger than the isothermal case shown 
in Fig. 1(b). Figure 3(b) has the parameter a decreased to 1 
m"' spreading the region of significant temperature van* 
ation'upward and decreasing the temperature gradient. The 
parameter /I is held constant at ten but since a was decreased 
this corresponds to a decreased frequency as is apparent 
from the interference pattern. The turning point for this case 
is at a height of about 1.6 m. 

IV. CONCLUSIONS 

The situation investigated here, plane waves approach* 
ing a finite impedance boundary at an arbitrary angle 
through a lapse temperature gradient is seen to demonstrate 
the major phenomena occurring at large distance from a 
point source. The interference pattern is seen to be signifi* 
cantly distorted from the isothermal situation by an inhomo* 
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geneoua layer about O.S m thick (a « 2 m~') and with a 
temperature increase of about 6*C near the ground (AT/ 

r-0.02) 
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